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We develop a general optimization strategy for performing a chosen unitary or non-unitary task 
on an open quantum system. The goal is to design a controlled time-dependent system Hamiltonian 
by variationally minimizing or maximizing a chosen function of the system state, which quantifies 
the task success (score), such as fidelity, purity, or entanglement. If the time-dependence of the 
system Hamiltonian is fast enough to be comparable or shorter than the response-time of the bath, 
then the resulting non-Markovian dynamics is shown to optimize the chosen task score to second 
order in the coupling to the bath. This strategy can protect a desired unitary system evolution from 
bath-induced decoherence, but ca also take advantage of the system-bath coupling so as to realize 
a desired non-unitary effect on the system. 
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I. INTRODUCTION 

Due to the ongoing trends of device miniaturization, 
increasing demands on speed and security of data pro- 
cessing, along with requirements on measurement pre- 
cision in fundamental research, quantum phenomena are 
expected to play an increasing role in future technologies. 
Special attention must hence be paid to omnipresent de- 
coherence effects. These may have different physical ori- 
gins, such as coupling of the system to an external en- 
vironment (bath) [H or to internal degrees of freedom of 
a structured particle , noise in the classical fields con- 
trolling the system, or population leakage out of a rel- 
evant system subspace [Sj. Formally, their consequence 
is always a deviation of the quantum state evolution (er- 
ror) with respect to the expected unitary evolution if 
these effects are absent In operational tasks such as 
the preparation, transformation, transmission, and de- 
tection of quantum states, these environmental couplings 
effects are detrimental and must be suppressed by strate- 
gies known as dynamical decoupling [5|-l8|], or the more 
general dynamical control by modulation [9l4l3|. 

There are however tasks which cannot be implemented 
by unitary evolution, in particular those involving a 
change of the system's state entropy [13, [l^ . Such tasks 
necessitate a coupling to a bath and their efficient imple- 
mentation hence requires enhancement of this coupling. 
Examples are the use of measurements to cool (purify) a 
system [T6l-[l9| or manipulate its state [20l - |22i or harvest 
and convert energy from the environment |23l427| . 

A general task may also require state and energy trans- 
fer (28|, or entanglement 1291 of non- interacting parties 
via shared modes of the bath |30l . Isij which call for maxi- 
mizing the shared (two-partite) couplings with the bath, 
but suppressing the single-partite couplings. 

It is therefore desirable to have a general framework 



for optimizing the way a system interacts with its en- 
vironment to achieve a desired task. This optimization 
consists in adjusting a given "score" that quantifies the 
success of the task, such as the targeted fidelity, purity, 
entropy, entanglement, or energy by dynamical modifica- 
tion of the system-bath coupling spectrum on demand. 
The goal of this work is to develop such a framework. 

The remainder of the paper is organized as follows. In 
Sec. |TT] we state the problem formally and provide gen- 
eral expressions for the change of the score over a fixed 
time interval as operator and matrix spectral overlap. In 
Scc. lIIII we discuss a general solution in terms of an Eulcr- 
Lagrange optimization. In Sec. II VI we apply the approach 
to the protection of a general quantum gate, which re- 
quires minimizing any coupling to the bath, whereas in 
Sec.|V]we consider the complementary case of enhancing 
the system-bath coupling in order to modify the purity 
(entropy) of a qubit. Open problems are outlined and an 
outlook is presented in Sec. I VII Supplementary informa- 
tion is provided in Apps. IXlandlB] 



II. OVERLAP-INTEGRAL FORMALISM 

A. Fixed time approach 

Assume that a quantity of interest ("score") can be 
written as a real-valued function P{t) ~ P[g{t)] of the 
system state g{t) at a given time t. This might be, for ex- 
ample, a measure of performance of some input-output- 
device that is supposed to operate within a predefined cy- 
cle/gate time t. Depending on the physical problem and 
model chosen, extensions and generalizations are conceiv- 
able, such as a comparison of the outcome for different 
t (on a time scale set by a constraint) [2^, a time aver- 
age P = J dr/(T)P[g(r)] with some probability density 
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/(t) [32], or a maximum P{t) = maxT-g[o,t]-P[^(''')] |33| . 
Here we restrict ourselves to the "fixed-time" definition 
as given above. Our goal is to generate, by means of 
classical control fields applied to the system, a time de- 
pendence of the system Hamiltonian within the interval 
< T < t that adjusts P{t) to a desired value. In particu- 
lar, this can be an optimum (i.e. maximum or minimum) 
of the possible values of P. Assume that the initial sys- 
tem state ^(0) is given. The change to the "score" over 
time t, is then given by the first-order Taylor expansion 
in a chosen basis 
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where the expansion coefficients 



dP 



^ iP)r 



(1) 



(2) 



are the matrix elements (in the chosen basis) of a Her- 
mitian operator P, which is the gradient of P[g{t)] with 
respect to g at t ~ 0, i.e., we may formally write P = 
(VgP)^Q = {dP/dgjf^Q. In what follows, it is P which 
contains all information on the score variable. Note that 
the transposition applied in Eq. ([2]) simply allows to ex- 
press the sum over the Hadamard (i.e. entrywise) matrix 
product in Eq. ([T]) as a trace of the respective operator 
product PAg. 

Let us illustrate this in two examples. If P is the ex- 
pectation value of an observable (i.e. Hermitian opera- 
tor) Q, so that P ~ Tr{gQ), then Eq. ^ just reduces 
to this observable, P = Q. If P is the state purity, P = 
Tr(p2), then Eq. © becomes proportional to the state, P 
= 2g{0). Note that the score P is supposed to refiect the 
environment (bath) effects and not the internal system 
dynamics. 

Equation ([T]) implies that Ag and with it P are small. 
Hence Ag must refer to the interaction picture and a 
weak interaction, while P[g{t)] should not be affected by 
the internal dynamics [so that no separate time depen- 
dence emerges in Eq. ([T]), which is not included in the 
chain-rule derivative]. In the examples above, this is ob- 
vious for state purity, whereas an observable Q might be 
thought of co-evolving with the internal dynamics. 

Our starting point for what follows is simply the rela- 
tion P = Tr(PA^) with some Hermitian P, whose origin 
is not relevant. 



B. Role of averaged interaction energy 

Equation ([T]) expresses the score P as an overlap be- 
tween the gradient P and the change of system state Ag. 
In order to find expressions for P in terms of physically 
insightful quantities, we decompose the total Hamilto- 
nian into system, bath, and interaction parts. 



and consider the von Neumann equation of the total (sys- 
tem and environment) state in the interaction picture. 
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-Q^Qtotit) = -i[-f^i(i),ftot(0] 



(4) 



[Hi{t) = ul{t)H[^\t)UF{t), (S) denoting the Schrodinger 
picture and Uvit) = T+e-'/od^[^s(^)+^B]]^ j^g solution 
can be written as Dyson (state) expansion 

^tot(0 - ^tot(O) + (-i) / dti[Hi(ti), ^tot(O)] 

+ {-if [ dh r dt2[i?l(<l), [Hl{t2), etot(O)]] + . . . , (5) 
Jo Jo 

which can be obtained either by an iterated integration 
of or from its formal solution 

Qtot{t) = m)Qtomu}{t), ui{t) = T+c-'/odt'iJi(o^ 

(6) 

by applying the Magnus (operator) expansion 
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expanding in ([5]) the exponential Ui = Ui(Hcs) and sort- 
ing the terms according to their order in Hi. 

We assume that initially, the system is brought in con- 
tact with its environment (rather than being in equilib- 
rium with it), which corresponds to factorizing initial 
conditions ptot(O) = g(0) (8) Pb- The environment is in 
a steady state ge, [^B,^fB] = 0, so it is more adequate 
to speak of a "bath" . Tracing over the bath in Eq. ([S]) 
then gives the change of system state A^ = g{t) — g{0) 
over time t, which we must insert into Eq. ([1]). We fur- 
ther assume a vanishing bath expectation value of the 
interaction Hamiltonian, 



(i/i)B =TrB(^Bi?l) =6. 



(9) 



As a consequence, the "drift" term corresponding to the 
first order in Eq. ([5]) vanishes, and we only consider the 
second order term as the lowest non-vanishing order ap- 
proximation. Finally, we assume that the initial system 
state commutes with P, 



[e(0),P] =0. 



(10) 



Hit) = Hsit) + HB + Hi, 



(3) 



In the language of control theory, Tr[p(0)P] is a 
"kinematic critical point" [U if Eq. pUj) holds, since 
Tr[c'^^(0)e-'^P] = Tr[g(0)P] + iTr(i7[e(0), P]) + 0{H^) 
for a small arbitrary system Hamiltonian H. Since we 
consider g in the interaction picture, Eq. (|10p means that 
the score is insensitive (in first order) to a bath-induced 
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unitary evolution (i.e., a generalized Lamb shift) j30| . 
The purpose of this assumption is only to simplify the 
expressions, but it is not essential. Physically, one may 
think of a fast auxiliary unitary transformation that is 
applied initially in order to diagonalize the initial state in 
the eigenbasis of P. Modifications to be made if Eq. (fTOj) 
does not hold arc provided in App. [Bj 

To lowest (i.e., second) order we then evaluate Eq. ([!]) 
for the score change as 



P = t^[H,P]H), H^- drHiir), (11) 

t Jo 



where (•) = Tr[^tot(0)(-)]- This expresses the change of 
score in terms of the interaction Hamiltonian, averaged in 
the interaction picture over the time interval of interest. 
Our scheme is summarized in Fig. [TJ 
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FIG. 1: Our control scheme: a system is brought in contact 
with a bath over a fixed time interval < t < t during 
which the time dependence of the system Hamiltonian Hs{t) 
is chosen such that a given system variable P is adjusted to 
a desired value at the final time t. 



C. Spectral overlap 

Alternatively, Eq. pT|) can be written as an overlap 
of system-and bath matrices, which allows a more direct 
physical interpretation. To do so, we assume d dimen- 
sional Hilbert-space, and expand the interaction Hamil- 
tonian as a sum of products of system and bath operators. 
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in such a way that (Bj) =0, which ensures that Eq. (|9]) 
is satisfied (otherwise we may shift 13 j = Bj — {Bj)I, Hg 
= Hs + Considering Eq. in the interac- 

tion picture and expanding Sj(t) = X]fc ^jk{t)Sk in terms 
of [Hermitian, traceless, orthonormalized to Tr(S'jS'fe) = 
basis operators Sj, defines a (real orthogonal) ro- 
tation matrix e{t) in the system's Hilbert space, with 
elements 



(13) 



where (•)id = Tr[(i~^/(-)]. These elements of the matrix 
e(t) may be regarded as the dynamical correlation func- 
tions of the basis operators. Analogously, we define a 



bath correlation matrix $(t) with elements 
^,k{t) = {B,{t)Bk)^. 



(14) 



It contains the entire description of the bath behavior 
in our approximation. Finally, we define a Hermitian 
matrix T with elements 



Vk,^{[S,,P]Sk 



(15) 



where (•) =Tr[^(0)(-)]. The matrix F may be understood 
as a representation of the gradient P with respect to the 
chosen basis operators Sj. Finally, we define the bath 
and (finite-time) system spectra according to 
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dTe"^^e(T). 



(16) 
(17) 



This allows to express the score, Eq. (jlll) . as the matrix 
overlap 

P = jj^ dtidhTi[e'^{ti)^{h-t2)e{t2)T] (18) 

/OO 
dcjTr[ett(w)G(cj)et(w)r] (19) 
-OO 

/OO 
dujTr[Ft{uj)G{u)]. (20) 
-OO 



In Eq. ([20)) we have used the cyclic property of the trace 
to write the spectral overlap in a more compact form 
by combining the rotation matrix spectra et(a;) and the 
gradient representation F to a system spectral matrix 
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Ft(w) = -e((w)rej(w 



t 



(21) 



Analogously, Eq. (jlSp can be written in a more compact 
form by introducing the matrix 



i?(ti,i2) =e^(ii)*(ii-i2)e(t2). 



(22) 



Equation (|20p is as a generalization of Kofman and 
Kurizki [ssj and demonstrates that the change P over 
a given time t is determined by the spectral overlap be- 
tween system and bath dynamics, analogously to DCM 
[l3| . or the measurement induced quantum Zeno and 
anti-Zeno control of open systems [U [s^. The bath- 
spectral matrix G{lu) must be positive semi-definite for 
all oj. If the same holds for the matrix Ft{uj), then P 
is always positive. Below we will consider such a case 
where P reflects a gate error and the goal is then to min- 
imize this error. The spectral overlap Eq. (pn)) can be 
made as small as desired by a rapid modulation of the 
system, such that the entire weight of the system spec- 
trum is shifted beyond that of the bath, which is assumed 
to vanish for sufficiently high frequencies. Since this fast 
modulation may cause unbounded growth of the system 
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energy, a meaningful posing of the problem requires a 
constraint. 

In general, Ft{uj) is Hermitian but need not necessarily 
be positive semi-definite, depending on the choice of score 
as encoded in T. This reflects the fact that P can increase 
or decrease over t. Depending on the application, our 
goal can therefore also be to maximize P with positive 
and negative sign. In what follows we will consider the 
question how to find a system dynamics that optimizes 
the score. 



III. EULER-LAGRANGE OPTIMIZATION 
A. Role of control, score and constraint 

Our considerations in the previous section suggest to 
define our control problem in terms of a triple {f,P,E) 
consisting of a control /, a score P, and a constraint E. 

The control is a set of real parameters //, which have 
been combined to a vector /. These can either be tim- 
ings, amplitudes, and/or phases of a given number of dis- 
crete pulses, or describe a time-continuous modulation of 
the system. Here, we focus on time-dependent control, 
where the /;(t) parametrize the system Hamiltonian as 
= Hs[f{T)], or the unitary evolution operator U{t) 
= T+c-'Io f^s{r'} = U[f{T)]. A direct parametriza- 
tion of U avoids the need of time-ordered integration 
of its exponent. The U{t) thus obtained js^l can be 
then used to calculate the system Hamiltonian Hs (t) = 

Two explicit examples of the score P pertain to the 
fidelity of P with a given pure state Fqi = (^|^|^'), (for 
which P = 1^) or to the von Neumann entropy which 
we can approximate (for nearly pure states) by the linear 
entropy, 5* = — fcTr(^ln^) w 5l = k[l — Tr(^^)], [for which 
P = ~'2kg{0)]. The latter score can be used to maximize 
the fidelity with the maximally mixed state g I (for 
which 5'l becomes maximum), or to maximize the con- 
currence q.tA3) = ^2(1 - Trgl), = TrB|4'AB)(^'AB|, 
as a measure of entanglement of a pure state I^'ab) of a 
bipartite system. 

If a constraint is required to ensure the existence of a 
finite (physical) solution, its choice should depend on the 
most critical source of error. An example is the average 

speed with which the controls change, E = dr / (r) , 
which depend on the control bandwidth in the spec- 
tral domain. A parametrization-independent alterna- 
tive is the mean square of the modulation energy, E 
= /pdr ((Aff)^(T)).^, where (•)id refers to a maximally 
mixed state and hence to a state-independent norm, and 
AH is the difference between the modulated and unmod- 
ulated (natural) system Hamiltonians. 



B. A projected gradient search 

We want to find controls / that optimize a score P{f) 
subject to a constraint E{f). A numerical local opti- 
mization can be visualized in parameter space as shown 
in Fig. H 




FIG. 2: Left: Optimization of the score P{f) subject to a 
constraint E{f) in control space {/} by walking along the 
component SP± of the gradient SP orthogonal to SE. Right; 
Resulting iteration consisting in steps determined by a small 
parameter e which yields a local solution depending on the 
starting point /q. 

We start at some initial point /q for which i?(/o) is 
the desired value of the constraint. Simply following the 
gradient SP would maximize or minimize P, but also 
change E. To optimize P while keeping E constant, we 
therefore move along the projection of SP orthogonal to 
SE, i.e., along SP± =SP— ^^g^^f SE. Since the gradients 
depend on /, the iteration consists of small steps fn+i 
= /„ ± eSP±{fn), 1. Assuming that neither SP nor 
SE vanish, the iteration will come to a halt where SP± 
vanishes, because the gradients are parallel, 

SP = XSE. (23) 

This condition constitutes the Euler-Lagrange (EL) 
equation of the extremal problem, with the proportional- 
ity constant A being the Lagrange multiplier. Its concrete 
form depends on the choice of P and E. Since the solu- 
tions of the EL optimization represent local optima of the 
constrained P, we may repeat the search with randomly 
chosen /o a number of times and select the best solution. 
The gradients at each point /„ may be obtained either 
from a calculation based on prior knowledge of the bath 
or experimentally from data measured in real time. A 
discretization of the time interval < r < t then reduces 
the variational (5 to a finite-dimensional vector gradient 
V. 



IV. GATE PROTECTION WITH BOMEC 

A. Gate error as average fidelity decline 

A particular application of our formalism is decoher- 
ence protection of a given quantum operation by bath- 
optimal minimal-energy control (BOMEC) [l3,[33]. Con- 
sider the implementation of a predetermined quantum 



5 



gate, i.e., unitary operation within a given "gate time" 
t. It is sufficient to consider a pure input state j^f). In 
the interaction picture with respect to the desired gate 
operation and in the absence of bath effects, we should 
therefore observe at time t the initial state |^). The 
quantity of interest is here the fidelity {'^\g{t)\^) , and 
we use the projector P ~ g{0) ~ |^')(^'| as the gradient 
operator, so that Eq. (fTO|) is satisfied and Eq. (tTTj) gives 
the fidelity change as the score 

P = (^lA^I*) = -t^{{^\H^\^) - {^\H\^)^)^, (24) 

which is given by H defined in Eq. 

Since a quantum gate is supposed to act on an un- 
known input state, we need to get rid of the dependence 
on l^*). One possibility is to perform a uniform average 
over all I^P). We may apply 



, , , , , , Ti-AB + Ti-ATi-B , , 
(*|A|vI/)(*|i3|*) = -— — (25) 



d(d+l) 



1381 139I I which gives the average 

d 



P 



(26) 



where (•)id = Tr[d^^/ (g) ^b(-)]- In Eq. (^5)) we have used 
TrsiJ = 6, which corresponds to Tr^ = in Sec. HTCl 
Because of this and because of Eq. ([9]), (i?)B = (^^i)b = 
0, we have (H)^^ = 0, and Eq. (|26p also describes the 
variance Var(i/) = (H^).^ - {H)i^- On the other hand 

P = —2kg{0) that gives the change AS* of entropy S = 
— fcTr(plng) is (up to a proportionality factor of —2k) the 
same as the P used here to give the change of fidelity, we 
have AS'=: —2kP. If we define a gate error £ as the aver- 
age fidelity decline, £ = —P, with P given in Eq. (p6|) . we 
can summarize the following proportionalities: gate er- 
ror = average fidelity decline ~ average entropy increase 
(purity decline) AS" ~ square (variance) of the average 
interaction energy H: 

2k 



£ 



-P 



'id 



1 



Var(i7). 



(27) 

In the matrix representation of Sec. Ill CI the average over 
the initial states in the matrix T defined in Eq. ((TS|) , gives 
r = -^/ [using TT{SjSk) = d5jk and TrS^ =0], so that 

^^dTlJ ^dc.Tr[e,(c.)4(c.)G(a.)] (28) 

in agreement with [37| [except a different normalization 
Ti{SjSk) = 26jk leading there to a prefactor j^]- From 
the requirement that £ > must hold for any positive 
semi-definite matrix et(a;)e|(a;), we conclude that G(w) 
must be a positive semi-definite matrix for any lj. The 
task of BOMEC is then to find a system evolution U{t) 
(cf. the control examples in the previous section) that 
minimizes £, subject to the boundary condition that the 
final U (t) is the desired gate. 



B. Comparison of BOMEC with DD 

It is appropriate to compare the effect of dynamical 
decoupling (DD) [1, @ with that of BOMEC. DD does 
not change with the bath spectrum G{uj). With an in- 
creasing number of pulses, DD shifts the weight of the 
system spectrum F{uj) towards higher frequencies, until 
the overlap Eq. (PH)) has become sufficiently small. This 
is illustrated for two different numbers of pulses of pe- 
riodic DD (FDD, pulses periodic in time), in the upper 
row of Fig. [3] in the case of a ID single qubit modulation 
(i.e., all pulses are given by an arbitrary but fixed Fauli 
matrix) . 
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FIG. 3: System modulation spectra F{lj) generated by two 
methods of DD. Upper row: periodic dynamical decoupling 
with n TT-pulses (PDDn), lower row: Uhrig dynamical decou- 
pling with the same number of vr-pulses (UDDn), compared 
for n = 11 (left column) and n = 19 (right column) pulses. 

Aperiodic DD such as UDD suppress low-frequency 
components (to the left of the main peak) in the system 
spectrum, which retain the system-bath coupling even if 
the main peak of the system spectrum has been shifted 
beyond the bath cutoff frequency (Fig. [3]). The plots in- 
dicate that this suppression of low frequency components 
is achieved at the price of a smaller shift of the main peak, 
i.e., shifting the main peak beyond a given cutoff requires 
more pulses in UDD than in FDD. Note that optimized 
DD sequences with improved asymptotics exist [7[ , which 
we will not consider here. 

System modulation spectra obtained with BOMEC are 
shown in Fig. 3] 

The plot refers to a qubit subject to pure dephasing 
(i.e., Z-coupling) by a bath whose spectrum G{io) has a 
Lorentzian peak and low-frequency tail. The BOMEC 
optimizes U{t) simultaneously for 3D Fauli matrix cou- 
plings to the bath (Z, Y and X). The resulting system 
spectrum F{uj) is shown for different energy constraints 
E which are increased in small and equal steps. For low 
E, F{uj) has a single peak on the left of the bath peak. 
Increasing E causes a second peak of F{uj) to emerge on 
the right of the bath peak, which continues to grow, while 
the peak on the left diminishes, until for high E, only the 
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FIG. 4: BOMEC-minimization of the gate error for a sin- 
gle qubit TT-gate caused by pure dephasing with a given bath 
spectrum [G{uj), bold red line]. The Z-components of the ob- 
tained system modulation spectra Fiiui) are shown for energy 
constraints Ei — 0.1 + 4(i — 1), i — 1,2, ... , 101, (thin lines, 
blue to green) individually scaled to 1. 



right peak remains. Fig. 2] hence demonstrates that the 
spectrum F{uj) generated by BOMEC changes continu- 
ously as E increases, but avoids overlap with the maxima 
of G{uj) irrespective of E. BOMEC can therefore be su- 
perior to all forms of DD including UDD, especially if the 
bath has high cutoff but bandgaps at low frequencies. 
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FIG. 5: Evolution within < t < t for an optimized cool- 
ing with initial p = 0.25 of (a) effective and (b) instanta- 
neous controls (red, green, and blue graphs show x, y, and 
z component), (c) ground state overlap of the system state, 
and (d) linear entropy [green, red and blue graph show nu- 
merical integration of the Nakajima-Zwanzig equation, time- 
convolutionless equation, and second order approximation of 
the solution [J], which are indistinguishable in (c)] (system- 
bath coupling strength K=W~'^,uJo = ^,t=W,E^ 100) . 



V. PURITY CONTROL OF A QUBIT 

To give an example of the opposite case, where the 
goal is to maximize the system-bath coupling, we apply 
our approach of constrained optimization to the linear 
entropy S'l = 2[1 - Tr(g2)] of a qubit. [Note that here Sl 
has been normalized to 1 by setting the coefficient k = 
d/{d — 1) = 2, cf. Sec. IIIII ] We assume an initial mixture 

g(0)=p|l)(l| + (l-p)|0)(0| (29) 

of a ground (excited) state |0) (|1)), where < p < 0.5 
is related to S'l by p = (1 - ^/T^^S^)|2. With Sj = 
aj denoting for d = 2 the Pauli matrices, Eq. (P^)) can 

be written in terms of Hr, ~ ^(Ts as p(0) = ° I? = 

J,"^H^J„ -I- ^ |°^JflLr, , where B = — — !^ is the inverse tem- 
perature. Purity and temperature are hence related via 
the energy scale luq. Our goal is a constrained optimiza- 
tion of ASl, i.e., P= -4 ^(0) in Eq. (g]). Unlike the gate 
error Eq. (j27p . ASl can be negative or positive, which 
can be understood as cooling or heating, respectively. 

The time evolutions resulting from a minimization 
of ASl for the initial state Eq. ([2^ are illustrated in 
Fig. [SJ The fj shown in Fig. [5l[a) are defined by U{t) = 
e-iMT)a3Q--^f2{T}a2Q'-^fi{T)a3^ whereas the ojj shown in 

Fig. [51[b) are given by Hs{t) = X]j The chosen 

constraint E = ^ J^dTTi-{Hs — Ho)^{t) can be written in 

" t ■ 2 ■ 2 

terms of the fj a.s E=IJ^ dti [fi + f2 + (/s - ojo)'^ + 
2/i(/3 — Wo) COS /2]. The overlap between the evolving 



system state ^(t) (in the Schrodinger picture) and the 
ground state |0) shown in Fig.[5ljc) indicates the fast uni- 
tary system modulation through short time population 
inversions without significantly altering the state purity 
as verified in Fig. [HJd). This can be visualized as fast 
TT-rotations of the state inside the Bloch sphere, which, 
together with smaller rotations, here result in the final 
reduction of S-^lt) seen in Fig. [Sfd). Fig. [SJ^d) also con- 
firms that for the chosen time and coupling strength, dif- 
ferences between various methods of approximation are 
small. 

In contrast to gate protection, no initial-state averag- 
ing is performed here, i.e., Eq. is known. Conse- 
quently, as Fig. [S] shows, the relevant components Fj = 
{Ft)jj{u}) of the system modulation spectrum contribut- 
ing to the spectral overlap Eq. (PH)) depend on the ini- 
tial state ^(0) via the matrix F Eq. p^ . [We assume 
an uncorrelated bath, i.e., Gjkij-^) ~ for j and Gj 
= Gjj{uj).] This influence is clearly visible in case of a 
constant (unmodulated, i.e., free) Hamiltonian (middle 
column) , for which we set = ^ with the final t being 
in the order of the bath correlation time. Cooling (heat- 
ing) is achieved by realizing negative (positive) ASl via 
maximum negative (positive) spectral overlap, as shown 
in the left (right) column of Fig. [HI This is the opposite 
to system-bath decoupling, where the goal is to minimize 
the overlap. 

The plots illustrate the role of the energy constraint 
E: increasing E allows to establish overlap with higher 
frequency components of the bath spectrum. This also 
suggests that for a bath spectrum with a finite frequency 
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minimization (maximization) even for weak modulation 
(i.e., for small E), is to adapt the temperature of the 
bath such that for an undriven system Hamiltonian J^O) 
no change is observed, AS'l = 0, which is a necessary 
condition for a system-bath equilibrium. This is would 
require non- unitary system modulation, e.g., the effect of 
repeated measurements [40l - |42| . 



VI. SUMMARY AND OUTLOOK 



FIG. 6: Cooling and heating of a TLS by minimization (left: 
a, b, c) and maximization (right: g, h, i) of the change of linear 
entropy for a given bath spectrum (red, Gj) with j = 1, 2, 3 
denoting x, y and z component. The optimized system spec- 
tra {Fj , green) are shown for an energy constraint E = 10^ 
(left) and E — 10'^ (right) as contrasted to the unmodulated 
Hamiltonian Ho (i.e., for E = Q, middle: d,e,f) and different 
initial states with p = 0.001, 0.05, O.f, . . ., 0.45, 0.499. 



cutoff, increasing E beyond a certain saturation value 
will not lead to further improvement of the optimiza- 
tion, cf. the general considerations in App. [X] In the 
time domain, increasing E leads to more rapid changes 
in the physical Hamiltonian however, requiring higher 
resolution of the numerical treatment. On the contrary, 
for attempted cooling (heating) by minimization (maxi- 
mization) of AS'l, a given E may be too small to lead 
to negative (positive) AS'l. The obtained ASl may then 
be understood as "reduced heating" ( "reduced cooling" ) 
as compared to a ASl obtained with an unmodulated 
Hq. This is shown in Fig. [T] The figure also illustrates 



\ASl 




opt. heat. ^ 100 




opt. heat. — 1 




■ ..unmixi.,-E — 




^^pt^cootiT^EL^^^^^ 




opt.cool.,^; = 100 





0.05 0.1 0.15 0.2 0.25 0.3 0.35 



FIG. 7: Change of linear entropy in units of the system-bath 
coupling strength obtained by minimization (attempted cool- 
ing) and maximization (attempted heating) of AS'l under dif- 
ferent constraints E — 0,1, 100 as a function of the initial p 
for a bath as shown in Fig. [G] 

once more that the ^(O)-depcndence of the spectra shown 
in Fig. ini is accompanied by a e(0)-dependency of the 
achievable change ASl for a given bath. For a maxi- 
mally mixed state in particular {p = 0.5), the matrix F 
Eq. vanishes and with it ASl. 

A possibility to achieve negative (positive) ASl by its 



Peculiarity of the approach: In summary, we have con- 
sidered a way of finding a time-dependence of the system 
Hamiltonian over a fixed time interval such that a given 
system observable attains a desired value at the end of 
this interval. The peculiarity of our approach is that it 
relies on knowledge of the bath coupling spectrum and 
adapts the spectrum of the system modulation to it. This 
allows to adjust the modulation to bandgaps or peaks in 
the bath coupling spectrum. In contrast to dynamic de- 
coupling of system and bath, which can be achieved by 
shifting the entire system-modulation spectrum beyond 
some assumed bath cutoff frequency, an enhancement of 
the coupling requires more detailed knowledge on the 
peak positions of the bath spectrum. In this way, our 
approach may comprise suppression and enhancement of 
the system-bath coupling in a unified way for executing 
more general tasks than decoherence suppression. The 
same approach can also be applied to map out the bath 
spectrum by measuring the coherence decay rate for a 
narrow-band modulation centered at different frequen- 
cies 



43| 



As far as the controls are concerned, we here consider 
time-continuous modulation of the system Hamiltonian, 
which allows for vastly more freedom compared to control 
that is restricted to stroboscopic pulses as in DD 0-0] ■ 
We do not rely on rapidly changing control fields that are 
required to approximate stroboscopic 7r-pulses. These 
features allow efficient optimization under energy con- 
straint. On the other hand, the generation of a sequence 
of well-defined pulses may be preferable experimentally. 
We may choose the pulse timings and/or areas as contin- 
uous control-parameters and optimize them with respect 
to a given bath spectrum. Hence, our approach encom- 
passes both pulsed and continuous modulation as special 
cases. 

Open issues: An open issue of the approach is the 
inclusion of higher orders in the system-bath coupling, 
which becomes important for strong or resonant system- 
bath coupling, so that a pcrturbative expansion cannot 
be applied. This may be the case especially when this 
coupling is to be enhanced in order to achieve a non- 
unitary operation (e.g. cooling), since in this case an 
optimization of the coupling may take us out of the do- 
main of validity of the entire approach. 

Another concern regards the initial conditions. Here 
we have assumed a factorized initial state of the system 
and bath. This prevents us from taking into account 
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system-bath-intcractions that may have occurred prior to 
that time. In particular, if the system is in equilibrium 
with the bath, their states are entangled or correlated 

An immediate problem of both higher order coupling 
and system-bath correlations is that their consideration 
requires knowledge of the corresponding parameters. It 
may be difficult to obtain such data with sufficient ex- 
perimental precision. Moreover, its consideration renders 
the theory cumbersome and the intuition gained from the 
spectral overlap approach presented here is lost. A way 
out is offered by replacing the "open" iteration loop in 
Fig. [2] with a "closed" loop [11], where the calculation 
of the score, constraint and their gradients are based on 
actual measurements performed on the controlled system 
in real time rather than on prior model assumptions, i.e., 
knowledge of bath properties. Such closed loop control 
would allow efficient optimization, but at the cost of los- 
ing any insight into the physical mechanisms behind the 
result obtained. 

From a fundamental point of view, it is interesting to 
derive analytic bounds of a desired score (under a cho- 
sen constraint) and see if this bound can be achieved by 
means of some (global) optimization, i.e., if the bound 
is tight. The need for a constraint in such optimization 
is not obvious if the task requires coupling enhancement, 
especially when the bath spectrum has a single maximum 
(Appendix A). 
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Appendix A: Bound estimation 

Our goal is to give constraint-independent upper and 
lower bounds for the maximum change P = Tr(PAp) that 
can be achieved with a given bath and P under the condi- 
tion PU)) . We assume that e, F, and G{uj) are quadratic 
(d^ — l)-dimensional matrices. If max[TrG(w)] < oo, we 
can estimate P by using that Tr{AB) < Tr(A)Tr(B) 
for positive semi-definite matrices A, B, and applying 
Holder's inequality in the form of f^^du) \f{u})g{uj)\ < 



sup{|.g(w)|} /^dw|/(a;)|. Decomposing T = Ti - r2 
into positive semi-definite matrices and making use of 
J do; et'^ {uj)et{uj) ^ I we thus get 

- P2 < P < Pi, Pt^t sup [TrG(a;)] TrF,. (Al) 

This reveals that for given t and G{uj), the bounds Pi 
depend on g{0) and P via F. 

Appendix B: Non-commuting score 

If Eq. ((T0|) does not hold, the following modifications 
must be made. We denote by A± = {A ± A^)/2 the 
(skew) Hermitian part of a given matrix A = A+ + A 
for convenience. Equation must be replaced with 

P = 2Ref dti/ dt2([i?i(ii),P]i?i(t2)) (Bl) 
Jo Jo 

= P + fdtJ \\t2 Tr{[etot(0), P]i?i(t2)Hi(ii)}, 
Jo Jo 

(B2) 

where (•) = Tr[^tot(0)(-)] and P is defined in Eq. pT|). 
Equivalently, we can write 

P ^2fdtJ dt2Tr\R+iti,t2)T++R^{h,t2)T4 (B3) 
Jo Jo 

= P-2 / dti/ di2Tr[flt(ti,t2)F_], (B4) 
Jo Jo 

where P is given by Eq. together with Eq. (f^ . In 
the spectral domain, the analogous expression is 

/oo 
dujTr[Ft{uj)g{uj)] (B5) 
-00 

/OO 
dwTr[F_(w)g^(a;)], (B6) 
-00 

where P is given in Eq. ff20|) and 

POO 

GiLj)= / dte'"**(i), (B7) 
Jo 

is related to Eq. ^ by G{iu) = 2g+{uj). 
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